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Abstract 

We derive an a priori error estimate for the numerical solution obtained by time and space dis¬ 
cretization by the finite volume/finite element method of the barotropic Navier-Stokes equations. The 
numerical solution on a convenient polyhedral domain approximating a sufficiently smooth bounded 
domain is compared with an exact solution of the barotropic Navier-Stokes equations with a bounded 
density. The result is unconditional in the sense that there are no assumed bounds on the numeri¬ 
cal solution. It is obtained by the combination of discrete relative energy inequality derived in m 
and several recent results in the theory of compressible Navier-Stokes equations concerning blow up 
criterion established in m and weak strong uniqueness principle established in US- 
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1 Introduction 

We consider the compressible Navier-Stokes equations in the barotropic regime in a space-time cylinder 
Qt = (0, T ) x R, where T > 0 is arbitrarily large and R C A 3 is a bounded domain: 

d t 0 + div^pu) = 0, (1.1) 


d t (gu) + div x (gu <g> u) + V x p(g) = div x S(V x u), (1.2) 

In equations (ll.llll.2p g = g(t,x ) > 0 and u = u(f,x) E R 3 , t E [0, T), iESJ are unknown density and 
velocity fields, while § and p are viscous stress and pressure characterizing the fluid via the constitutive 
relations 


§(V x u) 


li l V iT u + V* u - — div x ul 


M > °, 


(1.3) 
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p € C 2 ( 0,oo) flC^O, oo), p(0) = 0, p'(g ) > 0 for all £ > 0, lim ^ ^ = Poo > 0, (1-4) 

£—>•00 gl 1 

where 7 > 1 . 

Assumption p\ 0) > 0 excludes constitutive laws behaving as g 1 as g —> 0 + . Error estimate stated in 
Theorem Id.81 however still holds in the case lim e ^g+ §r=r > 0 at the price of some additional difficulties, 
see m for more details. 

Equations fll.llll.2l) are completed with the no-slip boundary conditions 


_c_ 

CC 

3 

II 

0 

(1.5) 

and initial conditions 


£>( 0 , •) = go, u( 0 , •) = u 0 , go > 0 in 12 . 

( 1 . 6 ) 

We notice that under assumption (11.31). we may write 


diVa;S(V x u) = /uAu + ^-VzdiVzU. 

O 

(1.7) 


The results on error estimates for numerical schemes for the compressible Navier-Stokes equations are 
in the mathematical literature on short supply. We refer the reader to papers of Liu |24] , [25], Yovanovic 

eh, Gallouet et al. Cl¬ 
in [EJ the authors have developed a methodology of deriving unconditional error estimates for the 
numerical schemes to the compressible Navier-Stokes equations (II. 11 - 11.61) and applied it to the numerical 
scheme (13.5113.7p discretizing the system on polyhedral domains. They have obtained error estimates for 
the discrete solution with respect to a classical solution of the system on the same (polyhedral) domain. 
In spite of the fact that m provides the first and to the best of our knowledge so far the sole error 
estimate for discrete solutions of a finite volume/finite element approximation to a model of compressible 
fluids that does not need any assumed bounds on the numerical solution itself, it has two weak points: 
1) The existence of classical solutions on at least a short time interval to the compressible Navier-Stokes 
equations is known for smooth C 3 domains (see Valli, Zajaczkowski [2Hj or Cho, Choe, Kim [2]) but may 
not be in general true on the polyhedral domains. 2) The numerical solutions are compared with the 
classical exact solutions (as is usual in any previous existing mathematical literature). In this paper we 
address both points raised above and to a certain extent remove the limitations of the theory presented 
in [171 . 

More precisely, we generalize the result of Gallouet et al. m Theorem 3.1] in two directions: 

(1) The physical domain 12 filled by the fluid and the numerical domain 12 h, h > 0 approximating the 

physical domain do not need to coincide. 

(2) If the physical domain is sufficiently smooth (at least of class C 3 ) and the C 3 — initial data satisfy 

natural compatibility conditions, we are able to obtain the unconditional error estimates with 
respect to any weak exact solution with bounded density. 

As in a, and in contrast with any other error estimate literature dealing with finite volume or mixed 
finite volume/finite element methods for compressible fluids (Yovanovich [31], Cances et al [5], Eyrnard 
et al. [D], Villa, Villedieu EDI, Rohde, Yovanovich m, Gastaldo et al. m and others) this result does 
not require any assumed bounds on the discrete solution: the sole bounds needed for the result are 
those provided by the numerical scheme. Moreover, in contrast with m and with all above mentioned 
papers, the exact solution is solely weak solution with bounded density. This seemingly weak hypothesis 
is compensated by the regularity and compatibility conditions imposed on initial data that make possible 
a (sophisticated) bootstrapping argument showing that weak solutions with bounded density are in fact 
strong solutions in the class investigated in HZ]. 

These results are achieved by using the following tools: 
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(1) The technique introduced in [T?J modified in order to accommodate non-zero velocity of the exact 

sample solution on the boundary of the numerical domain. 

(2) Three fundamental recent results from the theory of compressible Navier-Stokes equations, namely 

• Local in time existence of strong solutions in class (12.1111X171) by Cho, Choe, Kim [2]. 

• Weak strong uniqueness principle proved in m (see also HU)- 

• Blow ud criterion for strone solutions in the class (I2.11H2.12I) bv Sun. Wane. Zhane m- 

The three above mentioned items allow to show that the weak solution with bounded density 
emanating from the sufficiently smooth initial data is in fact a strong solution defined on the large 
time interval [0,T). 

(3) Bootstrapping argument using recent results on maximal regularity for parabolic systems by Danchin 

[S], Denk, Pruess, Hieber [3] and Krylov [21'j. The last item allows to bootstrap the strong solution 
in the class Cho, Choe, Kim [X to the class needed for the error estimates in m provided a 
certain compatibility condition for the initial data is satisfied. 


2 Preliminaries 


2.1 Weak and strong solutions to the Navier-Stokes system 

We introduce the notion of the weak solution to system (ll.lHl.4p : 

Definition 2.1 (Weak solutions). Let go : Q —>• [0,+oo) and uo : 12 —>■ M 3 with finite energy Eq = 
/ n (i£>o|w 0 | 2 + H(eo)) dx and finite mass 0 < Mo = f Q go dan We shall say that the pair (g, u) is a weak 
solution to the problem (HUHUSD emanating from the initial data (go,Uo) if: 

(a) g € C wea k([0, T]; L a (Ll)), for a certain a > 1, g > 0 a.e. in (0, T), and u € L 2 (0, T; W 0 1,2 (fi; M 3 )). 

(b) the continuity equation ( 11 . 11 ) is satisfied in the following weak sense 


/ gipdx 
Jn o 



(j)dty>+ gu ■ V x tp) dxdt, Vr € [0, T], V</? € C^°([0,T] x Q). 


( 2 . 1 ) 


(c) gu e CweakdO, T]; L b (Ll- 1 M 3 )), for a certain b > 1, and the momentum equation (II.2p is satisfied in 
the weak sense, 

/ gu-ipdx\ = / (gu ■ dtp + gu ® u : \7ip + p(g) divyddxdt 

Jn lo Jo Jn y 1 

(fiSIu : V^^dxdt +(p + A)divudiv(/?j dxdt, Vr € [0,T], \/p € C“([0, T] X 12;M 3 ). 


/o Jn 


( 2 . 2 ) 


(d) The following energy inequality is satisfied 

'1 


J (^-g\u\ 2 + H(g)j dx +J j ^/i|Vu | 2 + (// + A)| divu| 2 ^) dxdt < 0 , for a.a. r G ( 0 ,T), (2.3) 


with H(g) = g 


<-ep(z) 


dz. 


h z* 


(2.4) 


Here and hereafter the symbol / gdx\$ is meant for / g(r,x) dx — / go{x)dx. 

Jvi Jn 
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In the above definition, we tacitly assume that all the integrals in the formulas czu-ra are defined 
and we recall that C W eak([0, T]; L a (Cl)) is the space of functions of L°°([0, T]; L a (Cl)) which are continuous 
as functions of time in the weak topology of the space L a (Cl). 

We notice that the function g > H(g) is a solution of the ordinary differential equation gH'(g) — 
H(g) = p(g) with the constant of integration fixed such that H( 1) = 0. 

Note that the existence of weak solutions emanating from the finite energy initial data is well-known 
on bounded Lipschitz domains provided 7 > 3/2, see Lions [23] for ‘large’ values of 7 , Feireisl and 
coauthors [T3] for 7 > 3/2. 

Proposition 2.1. Suppose the Cl C R 3 is a bounded domain of class C 3 . Let r, V be a weak solution to 
problem II.1111,61) in (0,T) x Cl, originating from the initial data 

ro € C 3 (12), ro > 0 in Cl, (2.5) 

\ 0 €C 3 (n-,R 3 ), ( 2 . 6 ) 

satisfying the compatibility conditions 

Volan = 0, V x p(r 0 )\ d n = div I ,§(V I Vo)|ar 2 , (2.7) 

and such that 

0 < r <r a.a. in (0, T ) x Cl. (2-8) 

Then r, V is a classical solution satisfying the bounds: 

ll 1 / r ’llc([0,T]xr2) + H r llc 1 ([o,T]xn) + ll^ v ^ r llc([0,T];L6(n ; i?3)) + \\dt,t r \\c([o,T]-L 6 (n)) < D , ( 2 - 9 ) 

ll V llc 1 ([ 0 ,T]xn;.R 3 ) + ll V llc([0,T];C 2 (H;.R3)) + l|dtVa/y||c([ 0 ,T];L 6 (fi;i? 3 x 3 )) + ||dt it V|| L 2 ( 0 ) T;Z, 6 (n)) < D, (2.10) 
where D depends on 11, T, r, and the initial data ro, Vo (via ||(ro, N^o ) an ^ m in 2 . e ^ ro(x)). 

Proof: 

The proof will be carried over in several steps. 

Step 1 

According to Cho, Choe, and Kim [2], problem (11.1111.61) admits a strong solution unique in the class 
r € C([0, T m ); W 1,6 (Cl)), d t r € C([0, T M ); L e (Cl)), l/r € L°°(Q T ), ( 2 . 11 ) 

V € C([0, Tm\\ W 2 ' 2 (Cl] R 3 )) n L 2 (0, T m \ W 2 ’ e (Cl-, R 3 )), d t V € L 2 (0,T m ; Wo 1 ’ 2 ^;^ 3 )). (2.12) 

defined on a time interval [0,7m) , where Tm > 0 is finite or infinite and depends on the initial data. 
Moreover, for any T^ < Tm, there is a constant c = c(Tf r ) such that 

IMlL«( 0 ,TJJ f ;W' 1 ' 6 (n)) + ll^ r ||L°°( 0 ,T^;L 6 (n)) + l|l/ r l \ l °°( q t ) (2.13) 

+ ll v llL°°(o,T^;W 2 ' 2 (a ; i? 3 )) + ll v llL 2 (o,r* 1 -;W 2 ' 6 (n;iJ 3 )) + II^V|| L 2( 0]T * /; ryi.2(Q)) 

< c (ll r o||iy 1 > 6 (n) + ll v o||w 2 . 2 (o))- 


Step 2 

By virtue of the weak-strong uniqueness result stated in m Theorem 4.1] (see also m Theorem 
4.6]), the weak solution r, V coincides on the time interval [0, Tm) with the strong solution, the existence 
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of which is claimed in the previous step. According to Sun, Wang, Zhang [26], Theorem 1.3], if Tm < oo 
then 

limsup ||r(f)|| L oo ( m = oo. 
t->T M - 

Since (12.81) holds, we infer that Tm = T. At this point we conclude that couple (r, V) possesses regularity 
(I2.11U2321) and that that the bound (12.131) holds with c dependent solely on T. 

Step 3 

Since the initial data enjoy the regularity and compatibility conditions stated in (12.5112.71) . a straight¬ 
forward bootstrap argument gives rise to better bounds, specifically, the solution belongs to the Valli- 
Zajaczkowski (see [2HJ Theorem 2.5]) class 

r G C([0,T];W 3 ’ 2 (H)), d t r G L 2 (0, T; W 2 ’ 2 (H)), (2.14) 


V G C([ 0, T\- W 3 ’ 2 (H)) n L 2 (0, T; W 4 ’ 2 (n ; R 3 )), d t V G L 2 (0, T; W 2 ’ 2 (Q; R 3 )), (2.15) 

where, similarly to the previous step, the norms depend only on the initial data, r, and T. 

Step 4 

We write equation (11.21) in the form 

dtV - ^div x §(V x V) = -V • V*V + ^V x p(r), (2.16) 

where, by virtue of (12.151) and a simple interpolation argument, V G C' 1+I '([0,T] x fl;i? 3x3 ), and, by 
the same token r G C l+U {[ 0, T] x Q) for some v > 0. Consequently, by means of the standard theory 
of parabolic equations, see for instance Ladyzhenskaya et al. [22], we may infer that r, V is a classical 
solution, 

dfW. V 2 V Holder continuous in [0, T] x O. (2-17) 

and, going back to dm 

dtr Holder continuous in [0, T] x Q. (2-18) 


Step 5 

We write 

S7 x dtr = —V X V • V x r — V • V 2 r — V x rdiv. x V — rV x div x V; 
whence, by virtue (12.141) . (12.171) . (12.181) . and the Sobolev embedding W 1 ' 2 ^ L 6 , 

dtr G C'([0, T]; VL 1 , 6 (H)). (2.19) 

Next, we differentiate (12.161) with respect to t. Denoting Z = c^V we therefore obtain 

d t z - idiv a .S(V*Z) + V • V X Z = dt div x §(V x V) - dtV ■ V X V + dt Qv^(r)) , (2.20) 

where, in view of (12.191) and the previously established estimates, the expression on the right-hand side is 
bounded in C([0, T]; L 6 (H; i? 3 )). Thus using the L p —maximal regularity (see Denk, Hieber, and Pruess 
0, Krylov [2TJ or Danchin [SJ Theorem 2.2] ), we deduce that 

d 2 t \ = d t Z G L 2 (0, T ; L 6 (H; R 3 )), d t \ = Z G C([0, T\; W 1 ’ 6 ^; R 3 )). (2.21) 

Finally, writing 

d 2 t r = —dtY ■ V x r — V • <9 t V x r — c^rdiv^V — r^div x V, 
and using (I2.19p . (|2.2ip . we obtain the desired conclusion 

d 2 t reC([ 0,T];L 6 (H)). 

□ 

Here and hereafter, we shall use notation a ~ b and a ~ b. the symbol a ~ b means that there exists 
c = c(H, T, /U, 7 ) >0 such that a < c 6 ; a « b means a ~ b and b ~ a. 
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2.2 Extension lemma 


Lemma 2.1. Under the hypotheses of Proposition I~P1 the functions r and V can be extended outside 
12 in such a way that: 

(1) The extended functions (still denoted by r and V ) are such that V is compactly supported in 
[0, T] x R 3 and r > r > 0. 

( 2 ) 

ll V llc 1 ([ 0 ,T]x_R 3 ;_R 3 ) + ll V [|c([ 0 ,T];C 2 (i<; 3 ;i? 3 )) + ll^^ 7 3 ; V llc([ 0 ,T];L 6 (R 3 ;R 3 x 3 )) + ll^t V llL 2 ( 0 ,T;L 6 (i? 3 )) (2.22) 

~ ll V Hcl([ 0 ,T]xn ; R 3 ) + H V llc([ 0 ,T];C 2 (n;iJ 3 )) + ll^ V ^V|| C ([ 0]T ] ;i 6(Q ;R 3x3)) + \\d^ t V|| L 2 ( 0 ,T;Z,6(fi))! 

(3) 

ll r llc' 1 ([0,T]xR 3 ) + II^V 3 : r|| C ([ 0 iT ] ;i 6 ( / j 3 ;H 3 )) + ||9 t 2 tr|| C ([o,T];L 6 (i? 3 )) (2.23) 

~ ll r llc 1 ([0,T]xn) + ll^ V x^||c([0,T];L6(n ; i? 3 )) + \\dt,t r \\c{[o,T]-,L<i(n)) + 

II V llc 1 ([0,T] xTi;R 3 ) + ll V llc([0,T];C 2 (n;i? 3 )) + ll^ V ^ V llc([ 0 ,T] ; L6(n ; R3x3)) + || d\ t V IU 2 (0,T;L6(O))! 

(4) 

dtr + diVa;(rV) = 0 in (0,T) x R 3 . (2.24) 

Proof: We first construct the extension of the vector field V. To this end, we follow the standard 
construction in the flat domain, see Adams [TJ Chapter 5, Theorem 5.22] and combine it with the 
standard procedure of ‘flattening’ of the boundary and the partition of unity technique, we get (| 2 . 22 |) 
Once this is done, we solve on the whole space the transport equation (I2.24|) . It is easy to show that the 
unique solution r of this equation possesses regularity and estimates stated in (12.231) . □ 

Remark 2.1. Here and hereafter, we denote Xr(R 3 ) a subset o/L 2 ((0,T) x R 3 ) of couples (r, V), r > 0 
with finite norm 

II (a V)|| Xt (r 3 ) = IMIcAQO,T]xi? 3 ) + II^V a; r|| (: ;([ 0 jT ] ;L 6 ( /? 3 . R 3 )) + ||9f i t^||c([0,T];L 6 (_R 3 )) (2.25) 

ll^llc 1 ([0,T] xi? 3 ;ii 3 ) + l|V||c([0,T];C 2 (i? 3 ;R 3 )) + 11^ V|| C ([ 0 iT ] ;L 6( /? 3 ;i? 3x3)) + 11 V11X, 2 (O^jZ, 6 (^3)) 

We notice that the first component of the couple belonging to A'^(R 3 ) is always strictly positive on 
[0, T] x R 3 and set 


0 < r = min (M)e[0T]xR 3r(t,a:), r = max (t3;)e[0T]xK 3r(f,,r) < oo (2.26) 

2.3 Physical domain, mesh approximation 

The physical space is represented by a bounded domain 12 C R 3 of class C 3 . The numerical domains 12/, 
are polyhedral domains, 

Pft, = U^gx-A', (2.27) 

where T is a set of tetrahedra which have the following property: If K n L 0, K L, then K n L is 
either a common face, or a common edge, or a common vertex. By £(K), we denote the set of the faces 
a of the element I\ E T- The set of all faces of the mesh is denoted by £\ the set of faces included in 
the boundary <912/, of 12 /, is denoted by £ ex t and the set of internal faces (i.e £ \ £ ex t) is denoted by £ m \■ 
Further, we ask 

V/, E <912/, a vertex => V/, E <912. (2.28) 
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Furthermore, we suppose that each K is a tetrahedron such that 

£[K\ ~ diam [K\ « h, (2.29) 

where £[if] is the radius of the largest ball contained in K. 

The properties of this mesh needed in the sequel are formulated in the following lemma, whose (easy) 
proof is left to the reader. 

Lemma 2.2. There exists a positive constant dn depending solely on the geometric properties of 317 
such that 

dist [x, 317] < dnh 2 , 

for any x £ 317/,,. Moreover, 

1(17/, \ 17) U (17 \ 17/, ) | ~ h 2 . 

We find important to emphasize that 17/, cjL 17, in general. 


2.4 Numerical spaces 

We denote by Q/,(17/,) the space of piecewise constant functions: 

Qh(fifc) = {«£ L 2 (n h ) I MK £ T, q\K G R}. 

For a function v in (7(17/,), we set 

vk = -rpj- J vdx for K £T and II^u(.t) = ^ vkIk(x), x € 17. 


(2.30) 


(2.31) 


A'eT 


Here and in what follows, Ik is the characteristic function of K. 

We define the Crouzeix-Raviart space with ‘zero traces’: 

V h , 0 (Sl h ) = {v£ L 2 ( 17 ft ), Vif G T, G Pi(if), (2.32) 

Vo- G S int , cr = if|L, f v\ K dS = f v\ L dS, \/a'£ £ ex t , [ ud5 = 0}, 

J CT J <J J cr 

and ‘with general traces’ 

14(14) = {v £ L 2 ( 17), \/K £ T, V\ K £ Pi (if), V<x G 4nt, a = K\L, [ v ]K dS = f tpdS}. (2.33) 

J cr J cr 

We denote by Hff the standard Crouzeix-Raviart projection, and n' ( 0 the Crouzeix-Raviart projection 
with ‘zero trace’, specifically, 


r lX ■ 4(17/,) -> 14(17/,), / U v h [<f>] dS, = I <t>ds x for all o £ £, 


TV 


n^o : (7(17/,) ->• 14(17/,), / n^ ; 0 [4] dS^ = I (j) dS^ for all a G 4nt, / n£ 0 [<£] dS^ = 0 whenever a £ £ ext . 


tV 


v 


If v £ W 1 ’ 1 (17/,), we set 


Vrr = 


vdS for a £ £. 


cr 


Each element v £ 14(17/,) can be written in the form 


v(x) = v a ip a (x), x £ 17/,, 

cr£S 


(2.34) 


(2.35) 
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(2.36) 


where the set { (p a }a&£ C 14(0/,) is the classical Crouzeix-Raviart basis determined by 

V(<7 ,<j') €S 2 , [ (p a dS = 5 a ,a'- 

Similarly, each element v E V), o(0/,) can be written in the form 

v(x) = ^2 Va<Pa(x), x E fl h . (2.37) 

0"€:£int 

We first recall in Lemmas 12.31 - 12.71 the standard properties of the projection n) ( . The collection of 
their proofs in the requested generality can be found in the Appendix of HZj with exception of Lemma 
12.81 and its Corollary 12.11 We refer to the monograph of Brezzi, Fortin [3], the Crouzeix’s and Raviart’s 
paper [B], Gallouet, Herbin, Latche US] for the original versions of some of these proofs. We present the 
proof of Lemma 12.81 dealing with the comparison of projections n) t and IlJ ( 0 that we did not find in the 
literature. 

Lemma 2.3. The following estimates hold true: 

ll n ft[0]IU<»(x:) + Iln^oMIlL-CG ~ lldll(2.38) 
for all K E T and f E C(K); 


I10- n h W\Wlp{k) ~ hS \N s <j}\\ L p {K . R d*), s = 1 , 2 , l<p<oo, (2.39) 

and 

11 V(0 — II^[i^])| < ch s 1 ll^ s< / > ll_LP(A' ; R dS ), s = l)2, 1 < p < oo, (2.40) 

for all K E T and (f E C S {K). 

Lemma 2.4. Let 1 < p < oo. Then 

J2 W\ h K\ p ~ M P LP{ n h y ( 2 - 41 ) 

( t££ 


with any v E 14(0/,). 

Lemma 2.5. The following Sobolev-type inequality holds true: 

IMli«(n h )~ E [ \^xv\ 2 dx, 

K£T JK 

with any v E 14,o(^4)- 
Lemma 2.6. There holds: 

/ q div [a] dx = q div v dec, 

Jn 

for all v E C' 1 (ri/ l ,M rf ) and all q E Qh(£lh)- 

Lemma 2.7 (Jumps over faces in the Crouzeix-Raviart space). For all v E V/i,o(0/,) there holds 

1 f r 


(2.42) 


(2.43) 


E 

ad£ 


h 


E / .|V,^| 2 dx, 


KgT‘ 


K 


(2.44) 


where [v\a,n a is a jump of v with respect to a normal n a to the face a, 
Vx E (7 — A | L E £int; [d cr,n a (®) — 


v\k(x) - v\ L (x) ifria = ria,K 
v\l(x) - v\k(x) ifria = ria,L 


( n a,K is the normal of a, that is outer w.r. to element K) and 

Vx E a E £ ex t, Hff,n ff (^) = v(x ), with n a an exterior normal to <90. 
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We will need to compare the projections 11^ and nj( 0 . Clearly they coincide on ‘interior’ elements 
meaning K £ T, KildLlh = 0. We have the following lemma for the tetrahedra with non void intersection 
with the boundary. 

Lemma 2.8. We have 


11 n)T [<^>] - n^ 0 [<^111 z.- ^ 11 (li^ [0] - n^ 0 [^>]) 11 z,- ^ su P 

0 -cA'ndn 

for any cj) € C(K). 


ifKeT, K ndn h ^<b, 
(2.45) 


Proof: We recall the Crouzeix-Raviart basis (12.361) and the fact that 11^ and 11^ 0 differ only in 
basis functions corresponding to a £ S ext . We have 


nH 0 ]-n 


v 

h,0 


’(*0 - 


£ 

cre£(K) nfext 


f)dS 


< c(K) 


L°°(K) 


sup 

cr££(K) nfext 


I L°° (a ); 


(2.46) 


and 


v*(n£M-n£ 0 M)|| LOO(K )<h 


< ch sup ||<^>|| L oo W 
aCKndfl h 


E Vx'Pa-n H>ds 

■, , S „ \G\ j (7 


cre£(K)r£ext 


E 

o-e£(A)nfext 


L°°(K) 


Fc 


L°°(K) 


The proof is completed by || Jf ae £(K)n£ ext Vx<P*\\l°°(k) < c(K)h 1 . □ 

In fact, in the derivation of the error estimates we will use the consequence of the above observations 
formulated in the following two corollaries. 

Corollary 2.1. Let cj) £ C 1 (ii 3 ) such that 4>\dQ = 0- Then we have, 

mlW - KM\\l~{K) = 0 if K <=T,K n d Qh = 0, (2.47) 


\K[<t>] -KoM\l«>(k) +h\\V x (nl[cl>}-nl 0 m\Lo°(K-,RZ) ~ ^HV^IUco^.flS), (2.48) 

if K eT h , Kn dn h + 0, dK £ 012. 

Proof: Relation (12.471) follows immediately from (|2.45p . as there is an empty sum on the right hand 
side for ‘interior’ elements (K D 912 ^ = 0). 

For any x £ <912/,, there exists y £ d 12 (and thus 4>(y) = 0) such that 

|</>(x)| < dist[z,y]||V^|| L oo (jR 3 ;jR 3) ~ h 2 \\W x cl)\\ L oo^ R 3. R 3), (2.49) 

where we used Lemma o for the latter inequality. The proof is completed by taking supremum over 
K £ Th and combining with (12.491) . Note that the mesh regularity property (12.291) supplies a uniform 
estimate of constants c{K) from the previous lemma, which enables to write the latter inequality in 

(Em □ 

Corollary 2.2. For any f £ C(R 3 ), 

llltfM - KoMlHK) ~ ^ 3/p ll^ll L cc(n h) > 1 < p < oo- (2-50) 
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Proof: Apply inverse estimates (see e.g. |201 Lemma 2.9]) to (12.451) . □ 

We will frequently use the Poincare, Sobolev and interpolation inequalities on tetrahedra reported 
in the following lemma. 

Lemma 2.9. 

(1) We have, 

11^ - V K \\ L P(K) ~ h\\Vv\\ L p( K) , (2.51) 



Vct G £{K), ||u - v a \\ LP{K) ~ h\\Vv\\ LP(K) , 

(2.52) 

for any v G W 1,P (K), 

where 1 < p < oo. 


(2) There holds 

K’ — v k\\lp*(k) ~ |Vv||lp(a:), 

(2.53) 


Vcr G £(A'), ||u - V ff || £*•(*■) ~ I|Vu||lp ( a-), 

(2.54) 

for any v G W 1,P (K), 

1 < p < d, where p* = . 


(3) We have, 

11^ — v k\\li(K) — c ^\\^ v \\lp{K-,m2)' 

(2.55) 


\\v~Va\\ L*(K) < c/l^HVu||£P(AT ; R d )> 

(2.56) 

for any v G W 1,P (K), 

1 < P < d, where j = | + 


We finish the section of preliminaries by recalling two algebraic inequalities 1) the ‘imbedding’ in- 

equality 

(Ew) 1/ '^(Ew) 1/ ' 

(2.57) 


i— 1 2=1 


for all a = (a \,..., aif) G M L , 1 < q < p < oo and the discrete Holder inequality 

X>IN < (Ekr) 1/9 (El^l P ) 1/P , (2.58) 

2=1 2=1 2=1 

for all a = (ai,..., a L ) G M L , h = (&i,..., b L ) G R L , | + | = 1. 

3 Main result 

Here and hereafter we systematically use the following abbreviated notation: 

4> = u^[4>], ct> h = n v h [ct>], 0 M = n£ o [0. (3.i) 

For a function v G C([0,T],L 1 (H)) we set 

v n (x) = v(t n ,x), (3.2) 

where to = 0 < t\ < ... < t n _i < t n < t n+ \ < .. . t/v = T is a partition of the interval [0, T]. Finally, for 
a function v G Vh{&h) we denote 

\7 h v(x) = E V x u(x)1a'(x), div /t v(x) = E div x v(x)lA'(x). (3.3) 

K&T KeT 
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In order to ensure the positivity of the approximate densities, we shall use an upwinding technique 
for the density in the mass equation. For q € Qh(^h) and u € ~Vh,o(^h\ ®i 3 )j the upwinding of q with 
respect to u is defined, for a = K\L e £- mt by 


up _ J QK if Ua ■ riaj< > 0 
Q(J 


(-IL if «cr • n aK < 0 


(3.4) 


and we denote 


Up k(q, u) = Y ( l l a >u <y • n a,K = Y (qk[uo- • n ^ K ] + + qd^o • n ^ K \ ), 

crGf (A")cnfint <t££(K) n£int 

where a + = max(a, 0), a~ = min(a, 0). 


3.1 Numerical scheme 

We consider a couple (g n ,u n ) = (g n <(At,h) ^ u n,(At,h)^ 0 j ^numerical) solutions of the following algebraic 
system (numerical scheme): 

Q n eQ h (n h ), g n >0, VL n eV hfi (n h ;R 3 ), n = 0,l,...,N, (3.5) 


n n _ 0 n ~ v 

Y \I<\ K ., K 4>K + Y ^2 \ a \Q*' up ( u * ■ n a,K)4> K = 0 for any € Qh(^h) and n = 1 ,... ,N, 

K&r K&Tae£{K) 

(3.6) 

E §(M - er'ar 1 ) • OK + E E ■ «**] ■ VK (3.7) 

K&T KeT <j££(K) 

- Y p(Qk) Y M*u • n<r,K + M Y /vu n :Vu dx 

KeT ae£{K) KgT Jk 

+ — Y / divu n divu dx = 0, for any v € 14 o(U; R 3 ) and n = 1,..., N. 

3 i7dr jK 

The numerical solutions depend on the size h of the space discretisation and on the time step At. For 
the sake of clarity and in order to simplify notation we will always systematically write in all formulas 
( g n ,u n ) instead of u n, ( At,h ^). 

Existence of a solution to problem (13.5113.71) is well known together with the fact that any solution 
(t? n )i<n<iv C (Qh{ty) N satisfies g n > 0 provided g° > 0 thanks to the upwind choice in (13.61) (see e.g. 
P3 ED]). 

Remark 3.1. Throughout the paper, is defined in where u is the numerical solution constructed 

in 

3.2 Error estimates 

The main result of this paper is announced in the following theorem: 

Theorem 3.1. Let U C i? 3 be a bounded domain of class C 3 and let the pressure satisfy \l.f | ) with 
7 > 3/2. Let {g n , u n }o< n <jv be a family of numerical solutions resulting from the scheme \3.5dS. 71) . 
Moreover, suppose there are initial data [ro, Vq] belonging to the regularity class specified in Proposition 
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2. 1\ and giving rise to a weak solution [r , V] to the initial-boundary value problem M.1W1.6 1) in (0, T) xLl 
satisfying 

0 < r(t,x ) < r a.a. in (0, T) x Cl. 

Then [r,V] is regular and there exists a positive number 

C = C (M 0 ,Eo,r, r, bVfcr] > II (frr, Vr, V , d t V , VV, V 2 V) , 

ll^i r llL 1 (0,T;Ll' , (n))> H^^ 7r llL 2 (0,T;i 67/5 ' 1 ' _6 (n;K 3 ))> II ^^UlL 2 (0,T;L 6 / 5 (fi;R 12 ))> 


such that 


sup / 
i<n<JVJonn h L 


^ n |u n - V(f n , -)| 2 + #(*?") - H\r(t n , -))(g n - r(t n , •)) - H(r(t n )) 


dx (3.8) 

+A t J2 [ |V h u n -V x V(t n ,-)| 2 dx 

l<n<N Jnnn h 

<c(VAt + h a + [ — V 0 | 2 + H(g°) — H '(ro)(g° — r 0 ) — i/(r 0 ))l dxV 

V J nnfn L 2 J / 


where 


2 7 - 3 ,3 1 , 

a = - if — < 7 < 2 , a = - otherwise. 

j 2 — / — ) 2 


(3.9) 


Note that for 7 = 3/2 Theorem 13.11 gives only uniform bounds on the difference of exact and numerical 
solution, not the convergence. 


4 Uniform estimates 

If we take <fi = 1 in formula (13.61) we get immediately the conservation of mass: 

Vn=l,...lV, f g n dx = f g°dx. 

J£lh JQh. 


(4.1) 


Next Lemma reports the standard energy estimates for the numerical scheme (13.5113.71) . see again 

mm- 

Lemma 4.1. Let (g n ,u n ) be a solution of the discrete problem \3.5\\3. 71) with the pressure p satisfying 
Then there exist 

8a G [min(p^,^),max(£i^,p2)], a = I<\L £ £ int , n = l,...,N, 
e n K l,n G [min(p^ _1 , p^),max(^ _1 , g £)], KeT, n = 1,... ,1V, 

such that 


y. 1^1 ( 2 ^l -1 ^ 2 + H(6k)) ~ l-^l(o^l' u ^| 2 + 


i<gT 


KeT 

m . .. 

+ At ^ ^ (^p |V x tt n | 2 dx + (/r + A) / |divu n | 2 dx^) 
n=lftT V ^ /lr 


Ik 


+ [o,Tif"'] + [»(£!?"'] + K^“'l + Ks!£ t| ] = 0. (4.2) 
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for all m = 1,..., N, where 


ufs-uf- 1 ' 2 


[D^t ul ] = £ £ Wk 1 ^ k -^ k 

n =1 K£T Z 


MT” 1 ] = E £ 5 " |2 > 

n= \ KeT 

P5L a e “'] = A t f; v wrfcS 

n=i o-=ii:|Lef int 
m 

[o^L^i] = a* £ £ kl»"(ffi) 


,n _ _n—1 12 
c 

2 

n _ ,,,77.^2 

2 

(& - £l) 2 




\ u a ■ n CT ,/<|. 


(4.3a) 

(4.3b) 

(4.3c) 

(4.3d) 


n=l a=K\Le£ in t 

We have the following corollary of Lemma 14.11 (see m Lemma 4.1, Lemma 4.2]): 
Corollary 4.1. Under assumptions of Lemma o we have: 

(1) There exists c = c{Mq,Eq) > 0 (independent of n, h and At) such that 


( 2 ) 


where 


£(q,u\z,v) = / (g\ 

Jn h v 


n — 1 & —i^|Z/E£i n t 



N r 

fc£ / |V*u’ 

n=l JK 

2 | 2 dx < c, 




(4.4) 


N 







k £ \\ U ' 1 L6(a h ;K 3 ) — C ’ 

n =1 




(4.5) 


su Pn=0,...Arlk”w n 

2 II <r 

c. 



(4.6) 


SUp n=0 ,...7vl|£' n ||L7(O fl ) < C 




(4.7) 

to the 

class f/iere 

is c = c(M 0 ,E 0 

L) 7 ") L»(Q t ;1 1; 

)) > o 

.,N, 






su Pn=0,...N£{Q n iU n \f 

(t n ),U(t n )) < c, 



(4.8) 

(q\u - 

v 2 + E(g\z)^jdx, 

E(g\z) = 

H(g)~ 


*)-#(*)■ 


,r, \p'\c\ r ,r]) > 0 suc h that 






-Qlf 

1 l5S>i} 

L [max{ ok-, £>l}] 2-7 


1 <■ 

TO.(T,k\ < c 

tfl € [1,2) 

, (4.9) 


At£ £ W\{Qk - Ql) 2 \K ■ n a, K \ < c */7>2 

Tt = 1 CT = K\L/^.Sint 
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5 Discrete relative energy inequality 


The starting point of our error analysis is the discrete relative energy inequality derived for the numerical 
scheme (13.5113.71) in [17} Theorem 5.1]. 

Lemma 5.1. Let ( g n ,u n ) be a solution of the discrete problem \3.5W3.1\) with the pressure p satisfying 
([Z 2 P- Then there holds for all m = 1,... ,N, 

y. 2 1-^1 (jSk\ u k ~ Uk\ 2 - Qk\ u k ~ U° K \ 2 ^) + ^2 \K\(^E(g^\r^) - E(g° K \r ° K )) 


k&t ' 


K&T 


6 


+ At J2 E (» l\V x (u n -U n )\ 2 dx+^ / \div(u n -U n )\ 2 dx) <^T,, 

n=l KeT ' w 6 K 


i =1 


Ik 3 Jk 

for any 0 < r n € Q h {Ll h ), U n € V hfi (Qh', ^ 3 ), n = 1 ,... ,1V, where 

m ,, II r 

Ti = At V y (n / V x L/ n : V x (t/ n - u n ) dx + “ / divt/ n div(t/ n - u n ) dx), 

m 

T 2 = £ \ K \Hk 


n= l ifeT 
m. 


n=lK&T ae£{K) 
a=K\L 


n=l KgT ae£(K) 
cr=I<\L 


n= 1 XeT 
m 


\K\, 

At 


n ttti—1 

K ~ U K 

( u n K - ] + U\ 

-+o 

<1 

V 2 

/TT n 

K’ up ( K 


'\p(Qk)[U™ 


Q n K ){H\r%) - H\r" K - 1 )), 

K up H\r n K ~ 

l )[<-n a , K }. 


— u 


K h 


l^ / cr,K\ i 


n—1 KgT cf££(K) 
a=K\L 


(5.1) 


(5.2) 


6 Approximate discrete relative energy inequality 

In this section, we transform the right hand side of the relative energy inequality m to a form that 
is more convenient for the comparison with the strong solution. This transformation is given in the 
following lemma. 

Lemma 6.1 (Approximate relative energy inequality). Let ( g n ,u n ) be a solution of the discrete problem 
\3.5d3. 71) . where the pressure satisfies |1./| ) with 7 > 3/2. Then there exists 

c = c(^M 0 ,E 0 ,r,r, \p'\ c i[r,r\, II (d t r, Vr,V,d t V,W) ||loo ( q t;R i 8 ), 

ll^ r llL 1 (0,T;LT' , (f2))i ll^^ r llL 2 (0,T;L 6 '?/57-6(n ; ]R 3 ))) > 0, 

such that for all m = 1 ,... , N, we have: 

j n (g m \u m -VZ 0 \ 2 + E(g m \r m ))dx-(g°\u° - V 0 hfi \ 2 + E(g°\r°))dx 

m . .. 6 (6.1) 

+ A|V a (« n -n,o)l 2 ds+5 / \dW(u n -Vl 0 )\ 2 dx)<J2 s i + K A t+G rn , 

n=l KeT V Jl< 3Jl< ’ f^i 
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for any couple ( r i v) belonging to the class (2.25\) . where 

m , a r 

Si = At£ £ (/W V x V n ht0 :V x (V n hfi -u n )dx+^ / divW^div^o - «") dx), 

Jk 3 J K ’ ' > 


n —1 


=a (5: x: ia'i*- 


n t/- 

-1 v ft.O.AT v h,0,K f 


n=1 KeT 
m 


At 


ty"h,0,K — U K ) I 


5 3=a^ e e 


\v\e, 


n, up 


n=l KgT a££(K) 


/~n,up 


_ ^n,up 


Vh, 0 ,a~V^ K )V n h Z-^,K, 


( 6 . 2 ) 


S 4 = -At E / P(^ n ) div dx > 


n= 1 


S 5 = At 2 / (f" - 


p'(f n ) 


re=l 


[^r] n dx, 


5 6 = -AtE/ 


ra=l 


a, r 


and 


\G m \ < cAtJ2£( e n ,u n r n ,V n ), \R% At \ < c{VAt + h a ), 


(6.3) 


71=1 


with the power a defined in \3. ( J\) and with the functional £ introduced in 

Proof: We take as test functions U" = VJJ 0 and r n = r n in the discrete relative energy inequality m- 
We keep the left hand side and the first term (term T\) at the right hand side as they stay. The trans¬ 
formation of the remaining terms at the right hand side (terms T 2 —Tg) is performed in the following steps: 


Step 1 : Term Ti- We have 


T 2 — T 2 ,l + -^ 2,1 + R 2 , 2 , with T 2 .\ — At E E \K\Q 1 K 

n= 1 A'eT 


n —1 

-1 v h,0,K ~ v h,0,K 


V 


- V 


At 


(n», 


K ~ U K 1 1 


(6.4) 


and 


where 


and 


R2,l = At E E R 2,l , ^2,2 = At E ^2,2, 

n=l A"eT n=l 


,A'_ \K\ n _i ( V h,o,K - V hfl, K ) 2 \K\ ,-rOEEZEVg^) 


p7i,n. _ 11 „n- 

^2,1 ~ 9 Qk 


At 


Qk 


At 


R n 2,2 = - E l*l(?£ 

A'eT 


- V 


71—1 


— 1 h,0,K h,0,K ( 71 — 1 71 


At 


( w A r ‘ - u k) ■ 


We may write by virtue of the first order Taylor formula applied to function f 1 —)• V(t,x), 




n —li 


Jh,0,AT 


At 


1 


1 


|/\| Jk LAt L J tn _ 


dfV(z, x)dz 


h, 0 


dx 


\K\ Jk LAt 


— / [$V(*) (®)dz dx < ||[^V] M || L oo (0)T;i oc (C;R 3 )) < ||$V 


^ t'n — 


-I /z-,0 


L°° (0,T;Z/°° (Q;R 3 )) 


where we have used the property (|2.38|) of the projection n) t 0 on the space Therefore, thanks 

to the mass conservation we get 


l^2,’i I A ’ |A|At||9iP|| Lcx) j 0)r . Loo ^ n . R 3)). 


(6.5) 
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To treat term R 22 we use the discrete Holder inequality and identity dUD in order to get 


1-^2, 2 I ^ At cMo\\dtV\\ 2 Lx ^ QT . wlt00 ^. S 3^ +cMq / 2 ( £ lAl^rVr 1 


— u 


n 12 

AI 


KeT 


1/2 


IIW 


L°°(0,T;L°°(Q;: 


whence, by virtue of estimate (14.21) for the upwind dissipation term (|4.6a|) . one obtains 

|A 2 , 2 | < V / A ic(M 0 ,E 0 , ||5tVll^oo(6.6) 

Step 2: Term T 3 . Employing the definition (16.41) of upwind quantities, we easily establish that 

T 3 = T 3) i + 1 , 

m m 

with T„ = At Y. E E kl# up (ur p -vS")-v;, 0 , K <-n„, K , R 3:1 = AtY: E 

n=l AeT a<E£(K) n=l ae£ int 

|T^n. _ T^n |2 |T/-n _ V n 

1 j ->ti,(7 1 1 n I h,O.K h,0,L\ r n i —. i i 77 , I /i,0,A /i,0,Al r n 1 + w r-u t ^ c 

and # 3)3 = \<j\q k -—— -—— [u a ■ n^ K } + + \a\g L -—— -—- [u a ■ n^ L \ + , Vex = K\L € £ in t- 

Writing 

n, 0) A - n,o, L = ino - nk+n,A - n+n - n- 

+n,* - vi +n - +in - n, 0 k * = k\l e £ int , 

and employing estimates (j2.47|) (if AT D <90^ = 0) , (|2.48|) (if K D <90/,, / 0) to evaluate the L 00 -norm of 
the first term, (12.511) then (12.401b — i and (12.521) after (12.401) i to evaluate the L°°-norm of the second 
and third terms, and performing the same tasks at the second line, we get 


consequently 

whence 


\\Vhfi,K ~ ^, 0 ,lIIl°°(AUL;R 3 ) — ch\\W | (ACUZ,;M 9 ) > 

RZ I < h 2 C ll W ll^( 0 ,T)xQ;R 9 )M(^ + Ql)\K\, Vu = K\L € £ in t, 

5/6 


(6.7) 


|A 3 ) i| < h c|| 


At I2{Y, J2 %i 

n =1 K^T <j€.£(K) 


i( X! X! + 

AeT o-=a|Lg£(A) 


X 


u 


n |6 


1/3 


1/2 


( 6 . 8 ) 


< h c(M 0 , Aq, ||VV|| l°°(q t; r 9 ))) 


provided 7 > 6/5, thanks to the discrete Holder inequality, the equivalence relation (12.291) . the equiva¬ 
lence of norms (12.411) and energy bounds listed in Corollary 14.11 

Clearly, for each face a = K\L G £j n t, u" • + it" • n a _i = 0; whence, finally 


Ai = At^E E M<S' up ( 

n=l Kg'T ct££(K) 


■u" ,up - V| 


1 , 0 , <r) ■ (n,0,A — n,0,<r) U CT ' n o-,A- 


(6.9) 


Before the next transformation of term T 3 ) i, we realize that 

V/Va - VJ* i0iff = [V// 0 - Vj*k + V£ A - - + VJ* - V^ iCT + [V£ - V£, 0 ] CT ; 

whence by virtue of ([2.47H2.451) . (j2.51H2.52D and (I2.40p c -i . similarly as in (16.71) . 

ll^k,0,A ~ ^,0,o-Hl°°(A;R 3 ) ^ c/t || Va, V || A oo (o,T ;L°°(Cl:R 3 )) ) a C A. 


( 6 . 10 ) 
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Let us now decompose the term T 33 as 


T37 — T3 2 + i?3,2, with i?3_2 — At ^ -R3/2, 


71—1 


<r| e r p 


(LX - *r p ) ' (vt A , - VI„_ k )K' v ■ n„,K, and 


n,2 = &Y.i: £ 

Ti—i KeT cte.S(K) 

R b=E £ M<c up (LX -6 ;'" p )-(^- v !,o,k)(“: 

K(zlT <j(E.£(K ) 

By virtue of discrete Holder’s inequality and estimate (16.lOjl . we get 


n _ ^7i,up 


¥!'(t.K • 


1-^3, 2 I A c||VV|| ioo(QT . R9) 


(E 12 h \ a \& 

KeT erg £(K) 


n, up 


C UP - V 


2 \ 1/2 


£ £ hkii*“ p r) 1/(27 ” ) ( £ £ h\o\ 

KeT <t(z£(K) KeT cr(z£(K ) 


71, up 

h.O.cr 


ul - u n »P 


9 \ 1/9 


where ^ + | = 1, 70 = min{ 7 ,2} and 7 > 3/2. For the sum in the last term of the above product, 


we have 


< c 


(£ £ 

KE.T a££(K) 


12 12 %i 

K&T erg £(K) 


n _ nil q 
° II Li(K$ 


■u”-^ up ‘< c y ]r /ti 

ifeT erg £{K) 


+ ^ ik-<" 9 


icgr 


L1(K -]f 




2'yn — 3 /_ 

<c^( £ ||V 


icgr 


? ,n||2 

IIl 2 (/c 


y/i 

;R 9 J ’ 


where we have used the definition (USD, the discrete Minkowski inequality, interpolation inequalities 
(l2.55fE3Ul) and the discrete ‘imbedding’ inequality (12.571) . Now we can go back to the estimate of R 
taking into account the upper bounds (14.41) . (j4.7H4.8p . in order to get 


n 

3,2 


1-^3, 2 1 A h a c(Mq,Eo, ||VV|| L oo ( q t;R 9)), 


( 6 . 11 ) 


provided 7 > 3/2, where a is given in (16.31) . 

Finally, we rewrite term T 37 as 

771 

13,2 = r 3>3 + i? 3 j 3 , with R 33 = At Y2 $3,3> 

71—1 

771 

r 3 , 3 = a t j2Z E m c up (vlZ -<’ up ) • (n 0 ,a - n,o, k)KZ ■ nand (6-12) 

71—1 KtzT a(z£(K) 

R h = £ £ M* up (LX - <C up ) • (n„,„ - v^,*) (<.r p - LX) ■ 

KeT a^£(K) 


whence 


1^3, 3 i <c(iiwii ioo(QTiR9) ) Atx; ^ >0 ). 


71—1 


Step 3 : Term T 4 . Integration by parts over each K € T gives 


T 4 = -AtJ2 E ip(ei<) diY xV^ 0 dx. 

n =1 KeT "~ 


/K 


We may write 


||div x (Vo,/i — Vh)\\L°°(K) < c/i||V x VlL°°( 0 ,T ; L°°(n ;R 9 )), 


(6.13) 


(6.14) 
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where we have used (I2.47fl2~l8l) . Therefore, employing identity (12.431) we obtain 


r 4 = T 4 ,1 + i? 4 ,1, T 4>1 = - A i E E p(8k) div*V n dx 


n=l Ker • 


K 


(6.15) 


i? 4jl = -At ]T ]T / p(eSr)div*(n,o - n) d *. 

n=l KeT JK 

Due to dm) and (14.71) . p(g n ) is bounded uniformly in L 00 (L 1 (D)); employing this fact and (16.141) we 
immediately get 

l-R4.1l < h c(E 0 , M 0 , ||VV|| it x=( 0:T;i o°(f2 ; ]R9))). (6.16) 


Step 4: Term T§. Using the Taylor formula, we get 


H'{r n K ) - = H"(r n K )(r n K - r Jr 1 ) - ^H"\r n K )(r n K - r^ 1 ) 


'M (•z&i \(~.n _n—1\2 


where € [min(r^ ,r)t),max(r^ , r)^)]; we bifer 

m ( r n \ r n _ n -1 m 

T 5 = T 5 , 1 + j R 5)1 , withT 5il = At^ £ 1-^1 ( r A~ — Sk) P K ,^5,i = At^ £ and 


n=l Ae7" 

1 (V™ — r" -1 ! 2 

i?5,f = ^|A-|i7 w (r^) ( * J (g£ - r£). 


’ A' 


n=l KeT 


Consequently, by the first order Taylor formula applied to function t i->- r(t, x ) on the interval (t n -i,t r 
and thanks to the mass conservation m 


l-R 5 . 1 l < At c(M 0 ,r,r,\i/\ c ifcfl,\\dtr\\ LO o(Q T) ). 


(6.17) 


Let us now decompose T^i as follows: 


m ,, // n \ 772 

To a = T 5j 2 + 7?5,2) with T52 = At ^ E! S r K - Pk) P r n K [d t r] n dx, -R 5 , 2 = At E and 


p n,A' _ 
-^5,2 “ 


n=l A'eT 

J(„n\ , „n „n—1 


n=l a:sT 




(6.18) 

In accordance with (13.21) . here and in the sequel, [^r] n (x) = dtr(t n , x). We write using twice the Taylor 
formula in the integral form and the Fubini theorem, 


I u n ,K\ _ 1 

|77 5 ,2 I - At 


'n / b n 


V (rR')^A-(^A- - r%) / / d t r(z)dzd.sdx 


< 


P'( r K) rtn 


r n K Jtn-l JK 


lex-rS-iaM*) 


< AU|| e ” - f”||„ (K ) /" ||^rW|| iy(I0 

' Jin— 1 


dxdzds 

d^ds. 


Therefore, by virtue of Corollary 14.11 we have estimate 

17^5,21 < At c(M 0 ,E 0 ,r,r , |p'| c i([rx]> II^HIz^o.TjLVin))- 


(6.19) 
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Step 5 : Term Tq. We decompose this term as follows: 


t 6 = r 6)1 + r 6 , i, i? 6 ,i = E with 

n=lKeT aeS(K) 

m 

T 6 , 1 = E E E kleSr^'^K" 1 ) - ^'(C -1 ))< • *V,AT, and 

n=lKeT a=K\Le£{K) 


R 


K _ 

1,1 — 


M (l?^’ up - e n K ) (R'(r n K X ) - ^(rr 1 ))^ • n a ,K, for u = K\L € £ int . 


We will now estimate the term R^ ,r [' K . We shall treat separately the cases 7 < 2 and 7 > 2. The ‘simple’ 
case 7 > 2 is left to the reader. The more complicated case 7 < 2 will be treated as follows: We first 
write 


«'r*i < ^iiva'(r)ii 1 » (Qrl „, ) kn ( ,r p - ^i[ [ max{e ^7 t} p- T )/ 2 + 


• n^x 


%(?>!} [ maX W ; £T.}] (2 7)/2 + l{^<l} V^^/W ■ n *, K \, 

where we have employed the first order Taylor formula applied to function x 1—)■ H'[r[t n -i,x). Conse¬ 
quently, the application of the discrete Holder and Young inequalities yield 


1^6, i| < Vhc\\S7H'(r) Hloc ( q t . R 3 ) x 

m 

At E ( E E k|/i[l{e->i}[max{^, 

n= 1 A'eT o-ef(A') 

1 


|«<r • ™a,A 


1/2 






( 5/ MM# up - ex? 

A'eT cr=A'|Ae£(A') 

< V / /ic||VH / (r)|| L oo ( -Q T;R 3)X 

m 5 

A‘E{[|n»l» + ( E. H'*(*)* <2 - 7) ) 5 ] ( E V M'jK-'VkI 6 ) 1 

KeT cre£(K) 


7,K\) 


1/2 


71=1 


KeT 


E E MM<?r p - *) 2 

A'eT CT=A'|Aef(AT) 




L[max{^, £> l }] 2 7 


zr + !{c"<i} 


“a • n <J,A: 


1/2 


< 


\//ic||VH , (r)|| Loo(QT;R 3 ) { [|H h |e + ( J2j a \ %at)® (2 7) ) *] ( EM^I 

T At 


u 


n 16 


n=l ~ A'eT 

n,up _ n\2" 


cr££ 


1/6 




E [ E E wiHez™ - ii2- 7 + i {^<nj \ u « ■ n ^\ 

n=i KeTa=K\Le£(K) [max{g/e, q l \\ 

<Vhc(M 0 ,E 0 ,r,r, \p'\c(\r,r]), llVrll^Q^RS)), 

where, in order to get the last line, we have used the estimate (14.911 of the numerical dissipation to 
evaluate the second term, and finally equivalence of norms (12.411b — « together with (14.511 and (14.71) . under 
assumption 7 > 12 / 11 , to evaluate the first term. 

Let us now decompose the term Tq, 1 as 

771 

r 6 ,i = r 6 , 2 + Re, 2 , with T 6j 2 = At J2 E E M ^K H "( r K~ 1 )( r K^ 1 ~ C _1 )[< • n a , K ], 

n=\ KeTa=K\Le£{K) 


% = a<EE E K:C. 

n=lKeJC cre£(K) 

>n,cr,K 


and 


R, 


6,2 


- TC- 1 ) - ^"(rJ-'KrJ- 1 - rr‘))[< ■ n„,„]. 


19 


























Therefore, by virtue of the second order Taylor formula applied to function H ', the Holder inequality, 
(I2.41|) . and (14.51) . (14.71) in Corollary 14.11 we have, provided 7 > 6/5, 

m 

1^6,21 < hc(\H"\ C fa 7] ) + |^ W |c'([r,r]))l|Vr|| LOO(QT . R 3) AtJ2 II^IUr ( n ?t )||w n || L 6 (n?i . ]R 3 ) 

n =1 

< h c(M 0 ,Eo,r,f,\p'\ c i i]ttn) , ||Vr|| Loo( Q T;R 3 ) ). (6.20) 


Let us now deal with the term T, 67 . Noting that J X7r n 1 dx = E M( r cr r K may 

write Tg 2 = T 6j3 + i? 6i3 , with 


a&£{K) 


T 6i3 = -A I,, Qk R ’\ r n K X )u n • Vr n_1 dx, 

n=\ KeT’ 


IK 


Re,3 = At £ E ie?cH"(rl- 1 )(u n - u n K ) • Vr n_1 dx 

n= 1 A'eT" r " 


IK 


+aEE E 

n=l KgT <t££(K) 


Consequently, by virtue of Holder’s inequality, interpolation inequality (12.551) (to estimate ||u n — u k\\ l Y 0 , k R 3 ^ 

by /j( 5 7o-6)/( 2 7°)||Y_ r u n || i 2 ( ft -. R 9 \, 70 = min{ 7 , 2 }) in the first term, and by the Taylor formula applied to 
function x e-)- r{t n - i,x), then Holder’s inequality and (l2.55fE3Ul) (to estimate ||u” — u aII L 7 o^ R 3 ^ ^y 

/l (570 -6)/(2 70 )|| VrU n|| L2 ^ ;R9))) we get 

1-^6, 3 1 < h b c(M 0 ,E 0 ,r,r, |p , |ci(fc,rl)l|Vr|| L ao ( Q T;R 3 ) ), b= 57 ° 6 , (6.21) 

provided 7 > 6/5, where we have used at the end the discrete imbedding and Holder inequalities (12.571 
12.581) and finally estimates (14.41) and (14.71) . 

Finally we write = X /4 + iigq, with 


m r rJ (r n 1 

T 6 ,4 = -At E E / ■ Vr n dx, 

n=lKeT JK r K 

m r. 

Re,4 = At x; E & (H"(r n K )Vr n - "(r^ 1 ) Vr"” 1 ) • u" dx, 

n =lKeT JK 

where by the same token as in (16.191) . 

I-^ 6 , 4 1 A At c(M 0 , Eo,r, r, |p IcHIrx])’ II Vr, ^ r llL°°(Q T; K4), II^Vrll^Q^.^/cs-y-e^jR 3 ))), 


( 6 . 22 ) 


(6.23) 


provided 7 > 6/5. 

We are now in position to conclude the proof of Lemma 16.11 we obtain the inequality (16.11) by 
gathering the principal terms (16.41) . (16.121) . (|6.15|) . (16.181) . (16.221) and the residual terms estimated in 
(16.51) . (16.61) . (16. 8[) . (16.111) . (16.131) . (16.171) . (16.191) . (16.201) . (16.211) . (16.231) at the right hand side Y^i=\ R of 
the discrete relative energy inequality m- 

□ 
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7 A discrete identity satisfied by the strong solution 


This section is devoted to the proof of a discrete identity satisfied by any strong solution of problem EO 
[T5D in the class ([2.91 - 12.l(Jj) extended eventually to M 3 according to Lemma 12.11 This identity is stated in 
Lemma 17.1 1 below. It will be used in combination with the approximate relative energy inequality stated 
in Lemma 16.11 to deduce the convenient form of the relative energy inequality verified by any function 
being a strong solution to the compressible Navier-Stokes system. This last step is performed in the 
next section. 

Lemma 7.1 (A discrete identity for strong solutions). Let ( g n ,u n ) be a solution of the discrete problem 
\S. -711.V.71) with the pressure satisfying | l.Jffj , where 7 > 3/2. There exists 

c = c(^M 0 ,E 0 ,r,r, \p'\ c i^, \\{d t r, Vr,V,d t V, W, VV)|| Loo(Qt . R 4 5) , 

ll < ^ 2r ll L 1 {0,T-,LT(p,))i II^EdlL 2 (0,T;L 6 r/57-6(n ; K 3 )), \\&t V, <9fVV _ || i 2( 0)T;i 6/5(n ; Ri2))) > 0, 
such that for all m = 1 ,... , N, we have: 


Si + T^h.At, — 0) 

i =1 


(7.1) 


where 


m r. r. 

si = At E E / v,n, 0 = v*(n, 0 -u n )d X +^ divn i0 div(n )0 - u n ) dx), 


n= 1 KeT 

m 


s 2 = Atj2J2 \ R \i 


n=i KeT 

m 


K At 


(n, 


* = Atyy e MrrfKZ 

n=lK&T ae£(K) 

m . 

= —At ^ j p{f n ) div V n dx, 

^5 = 0, 


u p 

a ^cr 


n ’ u P 1 . 1 V™ 


■{V^~V^ K )V n h Z-n a , K 


%=\ 


S 6 = -A tJ2 / y(r”)w n • Vr n dx, 


i=i 


and 

\'K At \<c(h^ + At), 

for any couple (r, V) belonging to 1 H. 25\) and satisfying the continuity equation 11.11) on (0, T) x M 3 and 
momentum equation 11.HI) with boundary conditions 11.51) on (0, T) x Ll in the classical sense. 

Before starting the proof we recall an auxiliary algebraic inequality whose straightforward proof is 
left to the reader, and introduce some notations. 

Lemma 7.2. Let p satisfies assumptions Let 0 < a < b < oo. Then there exists c = c{a,b) > 0 

such that for all g E [0, oo) and r E [a, b] there holds 


E{g\r) > c(a,b)(l R+ \ [a/2t2b ](g) + H 7 lR + \[a/ 2 , 2 b](H) + (q ~ rf 


[a/2,2b; 


where E{g\r) is defined in fl/.#|). 
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If we take in Lemma 1721 o = o n (x), r = r n (x), a = r, b = r (where r is a function belonging to class 
(I2.25|) and r, r are its lower and upper bounds, respectively), we obtain 


E{g n (x)\f n (x)) > c{r,r)(l R+ \ lt/2 ^{Q n {x))+{Q n y(x)l R+ \ l r / 2 , 2 r]{Q n {x))+(Q n (x)-'f n {x)) 2 l lt/2 , 2 r]{Q n {x)) 


(7 ' 2) 

Now, for fixed numbers r and r and fixed functions g n , n = 0 ,,N, we introduce the residual and 
essential subsets of 11 (relative to g n ) as follows: 


Kss = {xen -r< g n {x) < 2r}, N? es = Q\N ( 


n 

ess’ 


(7.3) 


and we set 

[g]ess{x) = g(x)l N n s (x), [g] Tes (x) = g(x)l N n s (x), x £ ft, g € L 1 (H). 
Integrating inequality UZ2D we deduce 




dx < £(g n , u n 


r n ,V n ), 


(7.4) 


for any pair (r, V) belonging to the class (12.251) and any g n € Qh{^h)i Q n > 0. 

We are now ready to proceed to the proof of Lemma 17.11 

Proof: Since (r,V) satisfies (11.11) on (0,T) x H and belongs to the class (12.251) . Equation (11.21) can 
be rewritten in the form 

rdt,V + rV ■ W -f ’Vp(r) — pAV — p/2 >V div V = 0 in (0, T) x H. 

From this fact, we deduce the identity 

5 

E^ = ^o, (7.5) 

i—1 

where 

m n 

n 0 = At y / (r n [d t V] n + rV n ■ VV n + Vp(r n ) - pAV n - div • (VJJ 0 - u n )dx, 
n fe \n v 3 ) 

m r. m r. 

T\ = —At y / (gAE n + ^VdivP”) • {V^ 0 - u n ) dx, T 2 = Aty r n [d t V] n ■ (V% fi - u n ) dx, 

n =1 “ h n =1 ' “ h 

7YI TYl 

T 3 = At,y [ r n V n • VV" • (Yl o - u n ) dx, % = At y f Vp(r n ) • V n hfi dx, 
n=1 Jtth n= l 

m 

r 5 = 0, T 6 = -Aty / Vp(r n ) • u n dx. 

n= l Ja h 

In the steps below, we deal with each of the terms TZq and %. 

Step 0: Term TZq . By the Holder inequality 

m 

\U 0 \ < |H h \H| 5/6 c(r, b1cfc,r],ll(5 t r-,Vr,V,VV,V 2 V)|| Loo{QT . ]R 43 ) At^(||w-|| L6(n , i) + ||n i ollL6(n h )) 

n= 1 


< h 5/3 c(M 0 ,E 0 ,r, 


|p / |c [r ,r],ll(^r,Vr,V,VV,V 2 V)|| ioo(0T . I 


where we have used (14.51) and Q2.47H2.4gj) , (I2.38[) . 


(7.6) 
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Step 1 : Term T \. Integrating by parts, we get: 


T\ — 71,i + TZ 1 , 1 , 

m ,, 

with Ti,i = a153 E / ■ V(n,o - ^ div n,o div(n >0 - u n )) dx, 

n=lI<eT JK 6 

and 77]^] = 7i + J 2 , with 

m « 

h=AtJ2J2 / f^V(F” - VJ* )0 ) : V(n,o - « n ) + ^ div(V” - VJ* >0 ) div(Vj* i0 - u")) dx, ( 7 . 7 ) 

n=l I<eT JK 6 


h = -At 53 E E / • VV" • (n,o ~ O + f div V"(VJ* >0 - «") • dS 

n=l A'gT a&£(K) Ja 6 

m „ 

= - At E E / K-vr 

n=\ a&S Ja 


T/ n — n n 
v h ,0 w 


+ ^divF r 

J CT,ncr O 


T/ n 

v h,0~ u 


n<r) d S, 


where in the last line n a is the unit normal to the face a and [-Jo- n<T is the jump over sigma (with respect 
to n CT ) defined in Lemma o 

To estimate I\, we use the Cauchy-Schwartz inequality, decompose V n — V^ 0 = V 71 — + VJJ — 

VJJ 0 and employ estimates (l2.40L .-o. (I2.47H2.481) to evaluate the norms involving V(V ri — VJJ 0 ), and 
decompose 0 = VJJ 0 —VJJ+VJJ use (l2.47H2~4gD . (I2.39 D .-i . (I4.4D . the Minkowski inequalityt to estimate 
the norms involving V( 0 — u n ). We get 

|7l| < h c(Mq, Eq, IIVV, V 2 '\7|| L oo(0,T;L°°(a;R 36 )))- 

Since the integral over any face a E £; n t of the jump of a function from Vh,o{^h) is zero, we may write 


h = MY 


n=1 


E 

T^int 


fin 0 


(w n - (VV r 


u n -V 


h, 0 



(div V n )r^j 




a,n<j 



d 5; 


whence by using the first order Taylor formula applied to functions x i->- VV n (x) to evaluate the 
differences VV n — ( \7V n ) a , div V n — [div V n ] a , and Holder’s inequality, 


|/ 2 | < At/ic||V 2 V|| ioo(QT . R 27 ) Y, E \/l a l^ 

n = l (Tgfint 


1 

Vh 


o, n _ T/ n 

U V h,0 


a,n cT 


m 

< At/ic||V 2 V|| LtX)(Qr;R 27 ) ^ 5 Z (l a \ h 


+ 


1 


71=1 crgfint 


u n -V n 


h, 0 


2 

L 2 (cr;]E 


L 2 (cr;R 3 ) 

)■ 


Therefore, 

|77i,i| < hc(M 0 ,E 0 , ||V,VV,V 2 V|| Loo( g TiR 39 ) ), 
where we have employed Lemma 12.71 (14. 4 p and ([2.47hET3%l) , (12.39D . 

Step 2 : Term 7i • Let us now decompose the term 7i as 


(7.8) 


72 = 72,1 + 7*2,1, 

m 

with 72,1 = A t ^ ^ 
ti= 1 A'gT 1 


n —1 


AT 


— y 

aT 


71—1 


(n,o - u n ) dx, 77 2j1 = a n 2 ,i > 


and 77”f = J (r n - r n ~ l )[d t \} n • (VJJ )0 - u n ) dx + j r n ~ l [\d t V] n - 


n=i A'gT 
yn _ yn- 1. 

At 


(n, 0 -w n )dx. 
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The remainder TZ^’i can be rewritten as follows 


q-j n,K _ 
'^2,1 — 


1 


„n— 1 


d t r(tr)dt [ dt V] n iVl 0 -u n )dx+— / r" 

J At Jk L Jt n -iJs 


(' tn r tn 


dt v (z, -)dzds -(yl fl -u n ) dx] 


JK L Jt n -1 
whence, by the Holder inequality, 

\Kf\ < At[(||r|| i 0 o (QT) + ||a^|Uo 0(Qr) )(||^T^||^ (QT;M 3 ) |^T| 5 / 6 ( ||^| U6( ^ ) + ||^ ) 0 || L 6 W ) 

+ ||^ 2 y ,1 || L 6/5(n ; M 3 ))(||w"j|L6( A ') + IIV^oll^e^)) . 

Consequently, by the same token as in (I6.19P or (16. 2311 . 

1^2,il <Atc(M 0 ,£?o,nlKa t r,V,5 ( V,VV)|| i0 o (QT . R i6 ) ,||a2V|| i2(0iT ;L 6 / 5 (r2;R 3 ))) ’ (^-9) 

where we have used the discrete Holder and Young inequalities, the estimates (12.381) . (!2.47fl!H4&j) and 
the energy bound (14.41) from Corollary 14.11 

Step 2a: Term 72 ,i- We decompose the term 72 ,i as 


72,1 = 72,2 + ^2,2, 


with 72,2 = At ^2 5Z 

n=l ifeT 


m /■ , v n - V n_1 

n—1 


’ K 


' K 


At 

y n v”- 1 


(n,o - ^ n ) dx, 77 2 , 2 = At ^ ^ 77"; 

n=i xeT 


n,X 

2 J 


/* I/_ V 

and = Jj r”- 1 - r^ 1 ) J • (VJ,„ - u") dx; 


therefore, 


1^2,21 ~ I 51 ^2,2 I — ^ C ll^ r llL 00 (Q T ;K 3 )ll^^llL°°(Q T ;]f 

KeT 


|,,ti — T/ n II 

“ V /i,ollz,6(P;IE 


Consequently, by virtue of formula ([4.511 for u n and estimates (12.381) . (|2.47fEO%p . 

|77 2 , 2 | < h c(M 0 , H 0 , ||(Vr, V , <9*V, W)|| L oo (Qt;R i 8) ). 

Step 2b: Term 72,2 • We decompose the term 72,2 as 
72,2 = 72,3 + 7^2,3, 


(7.10) 


■71—1 


with 7^,3 = At 5] 5Z 

n=l KgT' 


\r n — y 

1 V h,0,K v h,0,K 


K 


K 


At 


(n,o - «”) d®, 77 2 ,3 = At ^ ]T 77"f, 

n=l iv eT 


and A = 


+ 

+ 


/c-‘( 




y n _ yn-l_ 


n T/n-1 


v 


At 


y™ _ y”- 1 - 


At 


(n,o-^ n )dx 


At U 

yn _ y 71- 1 


y 71 _ y Ti—l. 


At 


/iJC 


(n,o-u")dx 


At J /i,k 


r y n _ v n_1 


At J h,0,K 


y(vi 0 ~u n )dx= i« + i« + i«. 


We calculate carefully 


\T k \ - 
1 “ At 


h r x 1 J K ( J ^ dtV ~ ^ dtV ( z ^ h ’°\ K ' ~ un ^ dz ) dx 
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S A 7« 


,n— 1 


'tn — l 


[d t V(z)] h - [dtV(z )] h ,o 


A 


I 


m ,0 _ ' U ' , 1 |lL 6 (A';K 3 )d 2 - 


Summing over polyhedra K eT we get simply by using the discrete Sobolev inequality 


E n < 

AeT 


—r n ~ 
A t K 


T {( e 


) 1/6 ( e ii[wz)k-[a,v( 2) ] M r r> 2 


AeT 


AeT 


6/5 


5/6 > 


- aT* 


n —1 


< 


'tn—l 

5/6 nt 


VI,0 - - [W*)ko 


L6/5(fi h ;M 3 ) 


l 6 / 5 (a ; r 3 f 

dz 


At 


in,0 ~ ' un |lL 6 (n h ;R 3 )ll^i^( 2: )llL° 0 (n h ;R 3 )d^ ? 


where we have used estimate (12.501) to obtain the last line. 
As far as the term itf is concerned, we write 


I T k I — 1 r n ~ l 

\h I - A /a 


r / r /*^7l -| r ptn -| \ 

/ ( / <W*)dz - / &V(z)dz ) • (u" - n,o)ds 

« VL 7/ n _l L 7/ n _l 


h 

-Xt T * 


n —1 


ttn — l 


V, 


$V(z) 


1,6/5 


l' u " _ 


where we have used the Fubini theorem, Holder’s inequality and (12.511) . (12.401) i . Further, employing 
the Sobolev inequality on the Crouzeix-Raviart space Vh(Qh) (12.421) . the Holder inequality and estimate 
(12.401) a= i, we get 


E 1^1 — /T7 r A — ^fe,ollL6(O h ;R 3 ) 

AeT 


^71—1 


V^V^) 


A 6 / 5 (n h ;R 3 ) 


dz. 


We reserve the similar treatment to the term l [ Y . Resuming these calculations and summing over n 
from 1 to m we get by using Corollary 14.11 and estimates (I2.47H2.48D , (12.581) , 


1^2,31 A h °/ 6 c(Mo, -Em || (A V,W, dtV) || l°o(q t;R i6 ), ||3tVV|| l 2 (o,t ; l 6 / 5 (Q;R 9 )))- 
Step 2c: Term 72,3 - We rewrite this term in the form 


(7.11) 


77,3 — 77,4 + 7^2,4; 7^2,4 = At E E ^2,f> 

n=l AeT 


m-1 


with 77,4 = At ^ ^ 

n=l AeT 

V 


'A 


r „ rI v M j L _v MJi {u „ v „^ k)di _ 


and 7^’f = 


1 v h,0,K 


- V 


n —1 
h,0,K 


IK 


1 A 


At 


((«” - o - (n,» - di - 


(7.12) 


yn _y n_1 

First, we estimate the 7/ 00 norm of h ’°' K /± t h ’°' K as hr (10-51) . Next, we decompose 

n , 0 - no ,a = n,o - n+n - n,A+in - 

and use (l2.51L ,-o to estimate u n — u^, (12.511) .,,-^. (12.401) o—i to estimate V^' — R ) 3 K and (J257HZ5BD to 
evaluate ||[V)J — ^ 
we finally deduce 


evaluate ||[VJJ — ^/(,o]a|Ii / oo(A'R 3 ) — II V\ ~ ^'fe,ollL°°(A-R 3 )- Thanks to the Holder inequality and (14.41) 


172-2,41 < h c(M 0 ,E 0 ,r , \\(V, d t V , W)|| Z/00((3t . r i5 ) ). 


(7.13) 
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Step 3: Term T 3 . Let us first decompose 77$ as 
% = 73,1 + 7^3,1, 

m r. m 

with 773,1 = At E £ / r n K V n w • V V" • (VI o,K - O dx, ^ 3,1 = At £ £ , 

n= 1 K&T JK n =1 A'eT 

and = f (r n - r n K )V n ■ W n • (VJJ 0 - ix n ) dx + / r£(F n - V? „) ■ W n ■ (VJ* 0 - u n ) dx 

+ j K r n K (Vl 0 - n,o,A) • VV" • (n,o - «") dx 
+ J k r n K V n w ■ VV n • (v n hfi - Vl fitK - (u n - <-)) dx. 


We have 


by the Taylor formula, 


\r n - r^-|| L oo( A ') ~ h\\Vr n 


it 7-n \/-n n < 

I*' “ ^/i,0llz,°°(i^;M 3 ) ~ 


ivr 


by virtue of (12.391) ,-i and (|2.47U?7%8|) . 


It rn T/-n II ||T/-n T/n | 

I V /i.,0 V /i,0,AllL°°(A;R 3 ) — II v h,0 ~ V h, llL 00 (A;If 


Il°°(A;R 9 )> 


+ _ ^,aIIl°°(A;I 




H[v n h -vi 0 } K \\ LOO{K . R3) Zh\\vv r 

by virtue of (|2.51[l . (12.391) ,-i (12.4011 ,-1 and (j237HZ3SD ; 

ll w " — u aIIl°°(A;R 3 ) ~ ^II^ 7 ' u ' ,, ’IIl° 0 (A';R 9 )- 

Consequently by employing several times the Holder inequality (for integrals over K ) and the discrete 
Holder inequality (for the sums over K E T). and using estimate (14.41) . we arrive at 


1^3,i| < h c(M 0 ,E 0 ,r, ||(Vr, V, W)|| Loo(Qt . r i 5) ). 

Now we shall deal wit term 73 , 1 . Integrating by parts, we get: 


(7.14) 


/ w,o 1 rW' i -(n M -<)dx= £ w[v n h , 0 , K • na,*]vs ■ (n i0) A - o 

JK ct£S(K) 

= £ ki^[nw- n ^](K-nE0-(nv-<), 

o&£(K) 

thanks to the the fact that J2<tg£{K) fa Vh,K ' n o-,A'd5 = 0. 

Next we write 

m 

73,1 = 73,2 +^3,2, 7^3,2 = At E 7+3,2, 

n=l 


75,2 = a £ 

71=1 KdT aeS(K) 


W\C UP [VhZ ■ - n,A) • (VX - ^’ UP )’ 


(7.15) 


26 












and RJ 2 = E E MM - C up )[n,o,ir • n„, K ](v; - V’U) ■ (V" w - u%) 

KeT ae£(K) 

+ E E 


lair, 


n,up 


KeT ae£(K ) 


no,*: - C.) - n ’j< T" -*tu) - (v<u *o 


+ E E mc up [vX • n„, x ](K - n, K ) - («,« - w - w -*£;; 

KeT cre£(K) 


We may write 


K - ^h,o,K = V™ - V” + v” - v£ + v£ - v^ A + [V£ - V n hfi } K , 


and use several times the Taylor formula along with (l2.;->9lb -i . (12.511) . (I2.40lb -i . fl2.47H2.4Rl) fin order to 
estimate r A — f" ,up , V™ — 0 K , K — V^ p ) to get the bound 


1 ^ 3 , 2 I < hc\\r\\ Wl,oo(n)(l + || ^ |lu/ 1 .°°(Q T ;R 3 )) E '•l 


<7 \u 


KeT 


/ \ 2 

+c||r|| 

+ ll^lliyi,°°(Q T; R 3 ) J E E %H< 

KeT ae£{K) 


K\ 


- U„ 


We have by the Holder inequality 

E h W\\ u K\ < c ( E h \ 

KeT aeT 


fi a/6 

(7||W a | 6 1 < C 


El 

KeT 


\l, n — o, n II 6 
\ u u K\\lS(K$ 


1/6 


+(E 


KeT 


|U Il 6 (A';R 3 ) ) 


1/6 


< C 


(E « v “ 


KeT 


n \\L 2 (K-,R 9 ) 


1/2 


E E %IK - <1 < ■=[( E IK - <\\h ( K* 

KeT cre£(K) 


+(E E 

KeT ae£(K) 


i n _ n ||2 
** II r 9 r^.-irpS' 


KeT 
l/2-i 


1/2 


< /ic( E l|Vw„ 11 ^ 2 (^. 159 . 

KeT 


1/2 


viil 2 (a: ; r 3 ), 

where we have used (I 2 . 53 IL - 9 . (l2.51l - 1232]) r -9. Consequently, we may use (14. 4|) to conclude 

,W| 

Finally, we replace in 7i , 2 by Vl 0 a - V^ 0>A . We get 

m 

73,2 = 73,3 + ^ 3 , 3 , ^3,3 = At E ^-3,3> 


1^3, 2 I < hc(M 0 ,E 0 ,r , ||Vr, V, v v || £ oo ( Q T;R i5 ) y. 


n=l 


(7.16) 


m = At 2 £ x: ki^ up [C; 0 u E - n 0iA .). (vx - *r p ). ( 7 . 17 ) 

7 i=l KeT cre£(K) 


and 


^3,3 = E E M Tk - C UP )V n w • n CT , A ([F" - V M ]S - [VJI - n\ok) • (V«£ - <’ UP ), 

KeT ae£(K) 


committing error 


\ni 3 \ =< hc(M 0 ,E 0 ,r,\\Vr,V,VV\ 




)• 


(7.18) 
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as in the previous step. 

Step 4: Terms Ta We write 


r 4 = r 4 ,i + ft 4) i, 74,1 = - / Vp(r n ) • V n dx, 

Jflh 

7*4,1=/ Vp(r-).(V n -n,o)dx; 

Jft h 

whence 

1^4,l| < hc(r, |p , b|»; l r].||Vr|| £ oo ( Q T;R 3 ) ), 

by virtue of (j2.39|) —i . (J23ZH23SI)• 

Next, employing the integration by parts 


(7.19) 


7*4,2 - - 

KeT aE£ (K) ,<j€dQh 

Writing 


r 4 , 2 = t 4 , 2 + n A , 2 , r 4)2 = [ P {r n ) div v n dx, 

JQ h 

Jj(r n )V n • n^dS = - £ E jjj(r n )(v n - V^) • n^dS. 


Kd'T a££(K),(TEdQh ' 


v n - n i0i(7 =v* - n+n - n,.+m - v; 

we deduce by using &Ms=i, (1223.=!, (£52]) P =oo, J22ZD, (J23SD , 

ll-17-n T rn || < 

II V _ V /i,0,o-Hl°°(A';R 3 ) ~ 


n I 


Now, we employ the fact that 


whence 


£ £ /ds«i; 

KeTcrG£(K),a£dn h a 
|7*4, 2 | < Mg blcfcr], ||VV|| Loo(Qt; 


Finally, 


whence 


172-4,31 < Mb'lcfcr], ||(Vr, VV)|| Loo(Qt;R 12 ) ). 

Step 5 : Term 7f> We decompose T& as 

r 6 = r 6 ,i + 7 ^, 1 ,withr 6 ,i = - a t f; £ / p'(f> n • vnmx, 

n=\K&T JK 

m r. 

U^ = AtJ2J2 / (p'(r")-p'(r n ))-^-Vr"dx; 

Consequently, by the Taylor formula, Holder inequality and estimate (14.51) . 


(7.20) 


T 4 ,2 = T 4 ,3 + ^ 4 ,3, T 4i 3 = [ p{f n ) div V ri dx, 77 4 , 3 = f (p(r n ) - p(r n )) divV’Mx; (7.21) 

(7.22) 


(7.23) 


(7.24) 


|7*6,i| < /ic(Afo,fJo,r,r,|p , | C ri ( fc I r ]) >l|Vr|| LO o ( Q T;R 3 ) ). 

Gathering the formulae (17.71) , (17.121) , (17.171) , (17.211) , (17.231) and estimates for the residual terms (17.81) , 
(17.91 - 17. 13]) , (|7.14IT7TT81) , (17.191) , (17.201) , (17.221) , (17.241) concludes the proof of Lemma 17.11 □ 
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8 A Gronwall inequality 


In this Section we put together the relative energy inequality (16.11) and the identity (17.11) derived in the 
previous section. The final inequality resulting from this manipulation is formulated in the following 
lemma. 

Lemma 8.1. Let ( g n ,u n ) be a solution of the discrete problem 13.5W3. 7l) with the pressure satisfying 
where 7 > 3/2. Then there exists a positive number 

c= c(M 0 ,E 0 ,r,r, \p'\ c i\r,r], I \(d t r, Vr, V, d t V, W, V 2 V)|| I/ 00 (Qt;R 45 ) , 

\\dt r \\L 1 (0,T-,L'y'(a.))’ ll^^ 7 r ’llL 2 (0,T;L 6 r/5'y-6(Q ; R3^, \\d^V, dt W|| L 2( 0 ) T;L 6 / 5 (n ; R 12 ))) > 
such that for all m = 1,... , N, there holds: 


< c 


s(tr,« m \t m ,vZ>)+ A / iw(u”-v; 0 )i 2 dx 

Z n =lKeT JK 

m 

h a + VAt + S(e°, w°|r(0), V hfi (0))] + c At £ W, , Vh >0 ), 

71=1 


with any couple (r, V) belonging to 112.25]) and satisfying the continuity equation 11.11 ) on (0, T) x R 3 and 
momentum equation (ED with boundary conditions (ED on (0, T) x 11 In 1/ie classical sense, where a 
is defined in \6.3\) and £ is given in U-£\ )- 

Proof. We observe that 

772 r. 

Se - s 6 = At / p'(i 


r n — g n 


V" • Vr n dz +Al ]T / p\ 


72 = 1 ' 


71=1 


r n — n n 

r n ) —^-=_( u n — V”) • Vr n dx. 


Gathering the formulae (16.11) and (16.21) . one gets 


£{ e m ,u m \r m ,vl 0 )-£{g\u Q \m,V h M) + L^ £ |v(^-F^) 

n=1 K&T 

where 


I. 2 (A';l 


<£A+Q, (8.1) 


2=1 


72=1 KeT 
772 


72=1 KgT (j=K\Le£ K 

772 - 

v 3 = -At V / fp(e 


72=1 

772 


1 

s-cf 

n 

1—1 

1 

e 

?■ 

V h,Q,K (vn 

~ u k) j 

K J At 

' \ v h,0,K 

|a|(^ up -C’ UP ) 


• (no,. 

-p’(r n )(g n -r n )- 

-p(f n ))divV", 


r 11 - p n 

2 / n V n ' 

jin v , 

1 • \7r n dx, 



.€/ t 7-72 \ T/'^U.p 

A.,0,cr — v h,0,K) v h, 0,(7 ' n v,K ! 


Vi = Atj2 £ p'( f 

n=lK£T JK 

Q = ^hAt + Rh A t + G m . 

Now, we estimate conveniently the terms Vi, i = 1,..., 4 in four steps. 

V” _v n_1 

Step 1 : Term TT. We estimate the L°° norm of h ' u ' K ^ t h '°' K by L°° norm of dfV in the same 
manner as in (16.51) . According to Lemma 17.21 — r| 7 l R+ \ [ 2 -/ 2 ,2r ; ] C^) — c(p)E p (g |r), with any p > 1; in 

particular, 

\e - rf /5 l R+ \\rj 2 , 2 r](o) < cE(g\r) ( 8 . 2 ) 
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provided 7 > 6/5. 

We get by using the Holder inequality, 


E itfiw'-’V) 

KeT 


_ n —1 

—1 n—l\ v h,0,K v h,0,K 


V n — V 

* h n u * 


At 

1/2 


(VS. 


K ~ U K 


— c||9tW|| L oo(Q T;]R 3) 


3\ X 


( E WI &- 1 - ’•r 1 | 2 lfc/2,2f|(^)) + ( E - r 7T /5 lH + \fc/2,2F|((?A')) 


5/6 


A'eT 


KeT 


£l*l 

K&T 


V" _ n, n 

v h.O.A' u a: 


6\ 1/6 


<c(ll(ftV)|| I „ Wr;I1 . ) )(« 1 / 2 ( s”- 1 ,*”-i|i»- 1 ,<S” 0 I ) 

1/6 


+£ 5 /«(e”- i ,fi”- i if>"- i ,v;, 0 I ))( e im,o,A - -Aiii. (K;R .)) 


ATgT 


where we have used (18.21) and estimate (I4.8[) to obtain the last line. Now, we write V% 0 K ~ u \ = 
( \y"h 0 — u71 ]k — {V"h 0 — )) + (V/Jq — w n ) an d use the Minkowski inequality together with formulas 

IE33T1 . Vim to get ’ 


£ \\Vh,0,K ~ u k\ 


L 6 (K;R 3 ) 


1/6 


< (£ iiv(n l0 -« ; 


n\ i|2 


lL 2 (A";K 3 ) 


1/2 


K&T KeT 

Finally, employing Young’s inequality, and estimate OSD , we arrive at 

\V X \ < c(S,M 0 ,E 0 ,r,r,\\(V,VV,d t V)\\ Loo{QT ^ 5) ) 

m m 

A^ 0 ,w 0 |f°,v£)+At£^,n-r,v£))+Mt£ £ ||V(VJ* i 0 -w : 


n =1 


,n\ l|2 


n =1 XgT 


lL 2 (A';I 


(8.3) 


with any 5 > 0. 

Step 2: Term V 2 . We rewrite V^ a - V^ K = V^ - V\ K + [V^ 0 - + [V n hfi - VJ*]* 

and estimate the L°° norm of this expression by h|| W|| i 0 o/Q T . R 9 ) by virtue of (l2.47fEO%|) . f|2.51tf!Zln?l) . 
(I2.4 UI) o-i . Now we write V 2 = At )T// =] ^2 where Lemma 17.21 and the Holder inequality yield, similarly 
as in the previous step, 


[P 2 \ < c(r,r, ||VV’|| i 00 (Q T . R 9 ))x 

£ £ |ff|/i(f; 1 / 2 (ej up |r^ up ) + £ 2/3 (e? up |r^ up ) 


v n ’ up I I V n,up - fi^PI 

v h,0,cr \ I V h,0,cr “’cr 


KeT o-eE(K) 
< 


cfc?,||(v,vv)|| i . WTiR . I) )[( E E M'>( E « up IC up )) 

KeT cr££(K) 


1/2 


+ (£ £ W\hE(g^\C up 

Ka'T (t££(K) 


2/3-. 


(E E Mh\v n k Z-K Mr \) 

KeTaeE(K) 


6\ 1/6 


provided 7 > 3/2. Next, we observe that the contribution of the face a = K\L to the sums J^KeT 
Eae £{K) I a I hE(g™’ up |r” ,up ) and £a:gT T,cre£(K) W\ h \Vh^a-K^f is less or equal than 2|cr|/i(F;(^|f : ^)+ 
E(el\n)), an d than 2 |<j|/i(| V\ 0 K — u \| 6 + \V^ 0L — | 6 ), respectively. Consequently,we get by the 

same reasoning as in the previous step, under assumption 7 > 3/2, 


\v 2 \ <c(^,M 0 ,F;o,i,r,||(y,vv)|| L( » {QT;R i 2) )At£^,^r,F£)+Mt£ £ l|v(n, 0 -^)llV ;R 3 r 

n =1 AgT 

(8.4) 


n=l 
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Step 3 : Term V 3 . We realize that 


p(Q n K ) -p'(r n K )(g n K - r n K ) - p{r n K ) < c(r,r)E(g K \r K ), 
by virtue of Lemma 17.21 in combination with assumption (11.41) . Consequently, 

m 

\V 3 \ < c|| divV\\ L ~ {QT) At £ £{Q\u n \r\ V” 0 ). 

n=1 


(8.5) 


Step 4: Term V 4 . We write u n — V n as the sum ( u n — V^ 0 ) + (V ^ 0 — V n ) accordingly splitting V 4 
into two terms 


At E E L^'( 

n=lKeT 


r n — n n 
^n\‘ " 


r n — n n 
n\ 1 " 


IK 


(u n -V^ 0 )-Vr n dx and A P'( r ■ 

n =iKeT Jh r 


(V n K 0 -V n ).Vr n dx. 


Reasoning similarly as in Step 2, we get 

4 „/r m' 


|R 4 | < h c(S,M 0 ,E 0 ,r,r, |p'|c(fc,r])ll(Vr, VV r )|| i oo (n . I 

m m 

+ c(<5,||r,r,|p'| C ([^])l|Vr|| ioo(a;R 3 ) ) At £ £( e ", V" w ) + * At £ £ || V(V n hfi - u n ) \\ 2 L 2 {K . R 3 y 


n =1 


n=l ICeT 


( 8 . 6 ) 

Gathering the formulae m and (I8.3I) - (I8.6I) with 5 sufficiently small (with respect to g), we conclude 
the proof of Lemma 18.11 □ 


9 End of the proof of the error estimate (Theorem 13.11) 

Finally, Lemma 18. II in combination with the bound (14.81) yields 

m— 1 

£(e m , u m \r m , VZ 0 ) < c\h A + VAt + At + £(g°, u°\f (0), V hfi m] +cAtJ2 £(s n , u n \f n , vlo)- 


n=1 


whence by the discrete standard version of the Gronwall lemma one gets at the first step 

£ (f? m , u m \r m , V™,o) < c [ ha + VAt + £(e°, &°|f(0), V h , 0 (0)) 

Going with this information back to Lemma 18.11 one gets finally 

m r. 

£(g m ,u m \r m ,Vh, 0 )+At^Yl E / P V*(iC-V^ 0 )| 2 d X < c\h a +VAt+£(g° , u°|f(0), V\ 0 (0))1. (9.1) 


n=l KeT 


IK 


Now, we write 

Qk( u k ~ v h,o,A') 2 = Qk( u k - v?l ) 2 + ^QK^ n ( u K ~ v h,o,Ar) + Q n K (y n ~ v /P,o,at) 

where 

II V n — V' 


< 


< 


h,0,K\\L°°(K$ 

hi liv^w 


V n — v n 
l V V h IIl oo (A';]E 

^ + ||V X V?" 


+ livj*-y 


h,K\\L°°(K;l 


ll/°°(A';K 9 ) II v h\\L°°(K;S. 9 ) + ll^i — ^h,0 (A';R 3 ) ~ n ll V V 4 Il°°(A:;R 9 )' 

In the above calculation we have employed formula (12.391) to estimate the first term, estimates (12.511) o—i. 
(j2.40p .-i to estimate the second term, and formulas (I2.47P and (I2.48|) for K(ld£lh = 0 and K(idTlh / 0, 
respectively, to evaluate the last term. We conclude that 


+ \\m-V n h ,0\K\\L°°(K$ 

vv n| 


< 

3.3' ~ 


E l' U A ; _ V'k,0, K \ 2 - 8k\ u K - Vh,0,I<\ 2 ^) 


(9.2) 


ATeT' 
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>[ Q m (u m - V m ) 2 dx - [ q°(u° - V°) 2 dx + Li, 

Jnnfi h Jn nn h 

where 

\L\\ ~ h M 0 ||V ;c V|| L ^ ((0jT)xn;E 9 ) . 

Similarly, we find with help of (14.81) . 

\\ E (QK\r n ) -£(eSr,r n )|| L oo( K ) < hc{M 0 ,r,r, b| c i|r,?|l|Vr|| LO o { Q T;R 3 ) ); 

whence 

V \K\(E{ e n K \r n )-E{g 0 K \r 0 ) > f E{ 6 m \r m )dx - f E{g°\r°)dx + L 2 , (9.3) 

' JQnQh JQnQh 

where 

|L 2 | < h c(M 0 ,r,r, \p\ c i\r,r], l|Vr|| i oo ( g T . R 3 ) ). 

Finally, by virtue of (l2.47fEng]) and (12.401) s=2 

l|V(n,o - ^)IIl W 3) h\\(VV n ,V 2 V n )\\ Lao{K . Rl2 y, 

whence 

Aif E [ |V,(u--V^ 0 )| 2 dx> Atf; [ |(V/ l u n — Va;V n )| 2 dx + L 3 , (9.4) 

n=l KeT^ K n =l' Jnnn h 

where 

l-bsl < /t 2 c(||(W n ,V 2 V n )|| L oo (X;R 12 ) ). 

Theorem o is a direct consequence of estimate m and identities (19.2119.4|) . Theorem o is thus 
proved. 


10 Concluding remarks 

In the convergence proofs one usually needs to complete the numerical scheme by stabilizing terms, so 
that the new numerical scheme reads 


n n — n 71 - 1 

E \x\ 0K A t K <t>K+ E E w\qt v {K-^,k)^ k +tm) = ^ 


Ker h 


KGTh ad£(K) 

for any iG QhiS^h) and n = 1,..., N, 


( 10 . 1 ) 


E —( 

At v 


, Qk u k - Qk 1u k d ■ v k + E E 


XeT 


|a|^’U p <’ up [< • n*,K] • v* 


K^T cr(E:£(K) 


~ p{Qk) E kl^' n «T,if + M E / Vu n : Vu dx 


KgT 


ae£(K) 


KgT' 


K 


+— E / divu”divu dx + T m ((j)) = 0, for any u G 14 o(^( -R 3 ) and n = 1,..., N, 
3 ^T Jk 


where 


( 10 . 2 ) 


T c (</>) = h 1 ~ e E \°\le n U 

0 ‘C^int 


n a LV4I crn a » 


T m {4>)= E e € [ 0 , 1 ), 

0 ‘C^int 
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see Karlsen, Karper m, Gallouet, Gastaldo, Herbin, Latche m- These terms are designed to provide 
the supplementary positive term 

h 1 - £ 

to the left hand side of the discrete energy identity (321). They contribute to the right hand side of the 
discrete relative energy 32D by supplementary terms whose absolute value is bounded from above by 

/iM / 2 c(m 0 ,E 0 , sup ||r n ,U”,VU"|| Loo(n?i;R i 3 ) , sup sup [r n ] a ,n rT /h). 

' n=0,...,N n=0,...,N cr^Sint ' 

Consequently, they give rise to the contributions at the right hand side of the approximate relative energy 
inequality m whose bound is 


h (1 ' £)/2 c(m 0 , Eq, ||r, Vr, U, VU|| loo(Qt;R i 6) ) . 


Similar estimates are true, if we replace in the numerical scheme everywhere classical upwind formula 


33 ) 


Up K (q, u) = X q* Pu <T • n <r,K = X] (qxiuc r • n a , K } + + q L [u a ■ n ajK \ 


a&S{K) 

by the modified upwind suggested in m 


ae£(K) 

a=K\L 


Up K (g,u)= X ^-([ucT-n^K+h 1 £ ] + +[u a -n a: K-h 1 £ ] + ) +y {[u^-n^K+h 1 £ } +[u CT -n a , K -h l £ ] 


(10.3) 


o££(K) 
a=K\L 

where a = K\L € £; n t. We will finish by formulating the error estimate for the numerical problem 
fllO.II) . (110.21) or for (13.51) . (13.61) . (13.71) with modified upwind (110.31) . 

Theorem 10.1. Let U, p, [ro,V°], [r, V] satisfy assumptions of Theorem 1,5.11 Let (g n , u n ) n= o v .., 7 v be 
a family of numerical solutions to the scheme \3.5\) . Iil0.2\) or to the scheme 113. 5\) . \3. 6\) . \3. 71) 

with modified upwind UU.3\) . where e € [0,1). Then error estimate \3.8\) holds true with the exponent 


a = mm • 


•2 7 -3 1 — £ 
■ 7 ’~ 2 ~ 


}if§< 7 <2, a=^—^if 7 > 2. 
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